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Soliton-like solutions based on geometrically nonlinear Cosserat 

micropolar elasticity 

Christian G. Bohmer* and Patrizio NefP and Belgin Seymenoglu^ 


Abstract 

The Cosserat model generalises an elastic material taking into account the possible microstruc¬ 
ture of the elements of the material continuum. In particular, within the Cosserat model the 
structured material point is rigid and can only experience microrotation, which is also known 
as micropolar elasticity. We present the geometrically nonlinear theory taking into account all 
possible interaction terms between the elastic and microelastic structure. This is achieved by 
considering the irreducible pieces of the deformation gradient and of the dislocation curvature 
tensor. In addition we also consider the so-called Cosserat coupling term. In this setting we seek 
soliton type solutions assuming small elastic displacements, however, we allow the material points 
to experience full rotations which are not assumed to be small. By choosing a particular ansatz 
we are able to reduce the system of equations to a Sine-Cordon type equation which is known to 
have soliton solutions. 
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1 Introduction 

The theory of classical elasticity is based on the assumption that the material points are structureless 
and their possible internal properties like an orientation are disregarded. The Cosserat [4] model is 
a generalisation which takes into account a possible microstructure of the elements of the material 
continuum. In the Cosserat continuum material points can experience rotations without stretches. 
Therefore, in addition to the usual deformation field y there is an independent rotation field R. This 
idea motivated a rich variety of models with sometimes differing names ISHSmilHlllIJIMlIIIllIlIIl 
ElElllo]. 

Let us consider a material point and attach to it an arbitrary 3-frame. This means introducing 
nine new degrees of freedom which are 3 rotations, 1 volume expansion and 5 shear deformations. This 
general continuum is known as the micromorphic continuum [101125] . One can now formulate the special 
case where the structured material point is rigid and therefore it can only experience microrotation, 
this is the subject of Cosserat elasticity or micropolar elasticity which now has 6 degrees of freedom 
in total, three deformations ip and three microrotations R, considered in the present work. 
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This theory has two limits: In the case of no rotations one recovers standard elasticity theory, 
however, the limit of no deformations is also very interesting. Models of this type can in fact be traced 
back to MacCullagh, see [35], who was interested in understanding the ether. Since rotations in 3 
dimensions do not commute in general, the resulting theory of microrotations is generically always 
nonlinear. 

The Cosserat model has been theoretically proposed by the Cosserat brothers in a completely 
nonlinear setting [I]. However, they never proposed any specific constitutive relation. Their work was 
largely forgotten for nearly half a century when interest arose again in the late 1950s. Then, at first, 
scientists looked at a constraint model, in which the Cosserat microrotations are assumed to coincide 
with the continuum rotation. The investigations were mainly carried out in a linearized framework, 
see the recent acount given in [26] . Later, the full concept of free microrotations was rediscovered 
and soon extended along the given lines to micromorphic models, see [26] for a summary. It was 
quickly realized that these extended continuum mechanical models give rise to wave dispersion. Some 
intriguing new facets of linear micromorphic models being able to describe complete band gaps have 
only be discovered recently HlKIH]. Interestingly, this new feature really needs the micromorphic 
kinematics, stricly extending the Cosserat approach. The well-posedness for the static and dynamic 
case in the linearized setting is well-established. The static problem in a variational framework is 
also partially understood in the finite strain case, for both the Cosserat and the micromorphic model 
|2l|l7l|23]. 

The geometrically nonlinear micromorphic model was studied in [2211251128] where first existence 
theorems could be proved. To our knowledge, there are no known existence theorems for geometrically 
nonlinear Cosserat or micromorphic dynamic models. This point notwithstanding, in this paper we 
consider the geometrically nonlinear Cosserat model with a nonstandard constitutive coupling between 
Cosserat rotations and the deformation gradient. There are already many works taking nonlinear 
models with microstructure as a starting point for soliton-like solutions. However, most works seem 
to begin with a one-dimensional framework right away, while in our model we consider the fully 3- 
dimensional setting in the first place, see also [201 El [301 El [311 [38l [29]. It is only when we introduce a 
special ansatz that the Cosserat model simplifies. The simplification we use consists in assuming small 
displacements but allowing for arbitrarily large rotations. In doing so we keep the geometrically exact 
curvature expression and consider a certain nonlinear coupling between displacements and rotations. 

On the other hand, in exact wave-like solutions were found in the geometrically nonlinear 

Cosserat model. The existence of such solutions is a non-trivial fact, these solutions were hrst con¬ 
structed using spinor methods [1] which significantly simplified the Euler-Lagrange equations. This 
made it possible to find plane-wave solutions. Subsequent work [uia followed a more direct approach 
where wave-like solutions could be constructed without resorting to spinor methods. The study of wave 
type solutions in continua with or without micro-structure has always been of interest. The existence 
of linear elastic waves is well known and has motivated subsequent research ever since. 

A suitable deformation measure in Cosserat elasticity is the so-called dislocation curvature tensor 

_ 'Y' _ 

R Curli?, see m- It turns out that this commonly used naming convention is in fact slightly mislead¬ 
ing. When Cosserat elasticity is formulated in a differential geometric framework, then the orthogonal 
matrix R is interpreted as the frame or tetrad of a flat manifold with torsion. The object R diR is 
now a rank 3 tensor which is exactly the (Weitzenbdck) connection of this space, or the contortion 
tensor. The matrix R^CuilR is related to the Hodge dual of the contortion tensor, for more details 
see BE] and also [361137] . The axial part of this connection goes by various different names, wryness 
tensor [9], the Nye tensor m, and second Cosserat deformation tensor [10]. However, we will stick 
with the established notation of dislocation curvature tensor henceforth. 

In the above CurlM is the matrix curl, (CurlMj^^ = OkMuSkij, i.e. curl acting on the rows of 
M. Likewise we define DivM as the matrix divergence djMijei, again acting on rows. We will make 
frequent use of the Cartan-Lie decomposition of an arbitrary matrix in 3 dimensions 

AT = devsymM -I- skewM -I- - tr(M)l (1-1) 

o 


2 


where 1 is the identity matrix, symM = (M + M^)/2 is the symmetric part of M, skewM = 
(M — M^)/2 is the anti-symmetric part and dev M = M — l/3tr(M)l is the deviatoric or trace-free 
part of M. 

We also exploit the fact that every orthogonal matrix can be written as the matrix exponential of 
a skew-symmetric matrix, namely R = exp(A) where A is the skew-symmetric matrix generating the 
rotation. Furthermore, we follow the standard conventions of denoting by (p the deformation vector, 
u is the displacement vector and F = Vp = 1 -|- Vu being the deformation gradient. We define 
U = R F which is the non-symmetric right stretch tensor, also called the first Cosserat deformation 
tensor. In this context the second Cosserat deformation tensor is the above mentioned (Weitzenbock) 
connection R^diR. 

We remark that we write R and U with superposed bars in order to distinguish them from the 
factors R and U of the classical polar decomposition F = RU, m which R = polar(F") is orthogonal 
and U is positive definite, symmetric, which is the standard notation in elasticity. In particular this 
means U = VWf and R = polar(i^) = 


2 The complete dynamical Cosserat problem 


We write the energy functional for elastic deformation as 


Wiastic(F, R) = fi\\ symiR^F) - if + ^ tr(sym(i?^F) - 1)^, (2.1) 

where A and p are the standard Lame parameters. For the dynamical treatment, we will need to 
subtract kinetic energy which is of the form p/2||(pf. 

The energy functional of the micro-rotations VJ-otationai can be written in the form 

fcurvature(-R) = Ki || dev sym(i?^Curl i?)f -I- K 2 II skew(i?^Curl i?)f 


-I- K 3 tr(i?'^Curli?)^. 


( 2 . 2 ) 


Here Ki {i = 1,2,3) are three elastic constants. As before, in order to study the dynamical problem, 

we will need to subtract kinetic energy which is given by prot||.Rf = Pmt tr(i? R). 

Next, we wish to introduce a coupling between the elastic displacements and the micro-rotations. 
In order to do so, we will ‘couple’ the irreducible components of R Curl R with the corresponding 

—T 

irreducible parts of i? F — 1. This gives 


fdnteraction(-F, R) = Xi tr{R^Cur\R) tr(i?^F) 

-I-X 3 (devsymi?^Curli?, devsym(i?^F — 1)) (2.3) 


where xi and X 3 are the two coupling constants. In principle, we can also consider the skew-symmetric 

_ X _ j' _ _ Y _ T 

part of i? F — 1 in (skewF CurlF, skew(i? F — 1)). However, the term skew(F F — 1) is not used 
for the elastic energy (12.11) and therefore we will not consider it further. 

Finally, we will use the Cosserat couple term which is given by 

Hcoupiing(i^,^) = Pc llF^polar(F) - if, (2.4) 


where pc is the Cosserat couple modulus. When considering the Euler-Lagrange equations coming from 
this coupling term, we note that this will result in rather complicated expressions since we need the 
quantity 5polar(F)/5u. When considering small displacements, the situation simplifies signihcantly. 
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2.1 Small displacements and small rotations 

This nonlinear problem would be very difficult to address. Therefore we consider small displacements 

||u|| <C 1 and || Vu|| ^ 1 and small rotations || A|| <C 1 such that R = exp(A) = 1 + A + Here A is 

a skew-symmetric matrix, A = —A. In the following h.o.t. will stand for higher order terms. Then 

R^ F — 1 = exp(H)^(l + Vu) = (1 — A)(l + Vu) — 1 + h.o.t. 

= 1 — A + Vu — 1 -I- h.o.t. = Vu — A + h.o.t.. (2.5) 


Therefore we have that 


tr(i?^F — 1) = tr(Vu — A) -|- h.o.t. = tr Vu -|- h.o.t. 
sym(i?^F — 1) = sym(Vu — A) + h.o.t. = sym Vu -h h.o.t., (2.6) 

since tr(A) = 0 and also sym(A) = 0. Thus I4iastic in the linearised approximation is given by 

14iastic = /r||symVu|p -|- ^tr(symVu)^ (2.7) 

which is the classical elastic energy of linear elasticity. 

The term R Curl R simply becomes 

i?^Curli? = Curl A -I- h.o.t.. (2.8) 

Lastly, we need to consider polar(F). If = 1 + Vu, then one verifies that = F"’"F = 
l-fVu + Vu^-l- h.o.t. = 1-1-2 sym Vu -h h.o.t. and also that (1 -|- sym Vu)^ = 14-2 sym Vu 4 - h.o.t. 
so that we can write C/ = 1 4- sym Vu 4- h.o.t. which yields the result 

polar(F) = U~^F = (1 — sym Vu)(l 4 - Vu) 4- h.o.t. 

= (1 — sym Vu)(l 4- sym Vu 4- skew Vu) 4- h.o.t. 

= 1 4-skew Vu 4-h.o.t.. (2.9) 

Therefore we can write the Cosserat coupling term in the completely linear approximation as follows 

14oupiing = Me ||i?^polar(i^) - l|p = /ic 11(1 - A)(l 4- skew Vu) - l|p 4- h.o.t. 

= Me II1 — A 4- skew Vu — l|p 4- h.o.t. 

= Me II skew Vu — A||^ 4-h.o.t.. (2.10) 

2.2 A non-linear curvature model 

In the following we will consider a slightly different model which is different to the quadratic approx¬ 
imation as it is non-linear in the rotations. As before, we will assume that displacements are small 
||u|| <C 1 and ||Vu|| <C 1, but we allow for arbitrary large rotations R = exp A. However, we will 
assume that the terms VuA are small, this means ||VuA || <C 1 for all n > 1. 

In this case, we again find 


- 1 = Vu - A-4 h.o.t. (2.11) 

and our elastic energy is still ism . The curvature energy (1^ is unchanged and the interaction 
energy (12.31) is linearised in u, as before. The Cosserat coupling term requires a closer look. Namely 
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we find 


Me 

||i?'^polar(F) - if 


Me 

|i?'^(l -k skew Vu) — I p + h.o.t. 


Me 

|i?^ + R^ skew Vu — 1 p -k h.o.t. 


Me 

|i?^ -k (1 + A -k H^/2 -k .. .)'^ skew Vu — I p -k h.o.t. 

(2.12) 

Mo 

|i?^ — 1 + skew Vu p -k h.o.t. 


Me 

1 skew Vu — (i? — 1) p -k h.o.t.. 

(2.13) 


The crucial step is in (12.121) where we used the definition of the matrix exponential and our assumption 
that II Vu^ II ^ 1 but did not assume the rotations to be small. Combining this into a single energy 
functional gives us our non-linear curvature model 


V "2 


^elastic “t“ V^urvature “t“ Vinteraction T i^ouplingi 


(2.14) 


where the individual parts are given by 


V^eiastic = M II sym Vu||^ ^ tr(sym Vu)^ 

V'curvature = || dev sym(i?^Curl i?) || ^ -I- K 2 II skew(i?^Curli?)||^ 

-I- K 3 tr(i?^Curli?)^ 

^interaction ~ Xl tr(i? Curl i?) tl (Vu) 

+ X3 (devsym(i?^Curli?), devsym Vu) 

V'coupiing = Me II skew Vu - {R- 1)||^. 


The equations of motion are derived using the calculus of variations. We will vary the Lagrangian 
V 2 with respect to Sui and SAij. However, it is more convenient to work with the vector a = (oi, 02 , 03 ) 
such that Aij = EikjCLk and consider the variations with respect to this vector, we note SAij = SikjSak- 
We use the notation a = ymf + a| + a| for the norm of the vector a. The variations with respect to 
Sui are pretty straightforward to calculate. This cannot be said for the variations with respect to Sak- 
We note that the variations of the rotation matrix R with respect to Sak are already very complicated, 
the result is stated for convenience 


^Rij 

ScLk 


a cos a — sm a 


sin a 

ClkCll Silj ^ikj 


2 (cos a — 1 ) -I- a sin a 


a 

1 — cos a 


Hki^niCLj O.nQ’jiSij'^ 


{cliSjk Clj6ik ‘^Q'kSij^. 


(2.15) 


It becomes clear that the variations of V 2 with respect to Sak will result in very long expressions 
which we will not state explicitly. In particular, the variations of Wurvature will result in complicated 
expressions since this potential is fourth order in R. 

Let us now consider the variations of the functional V 2 with respect to the displacements Ui. This 
results in the following equation using the index notation 


pdttUi = (m + Me) Ui,jj + (m + ^ - Mo) Ujji 

+ ^(3Xi - X 3 ) Cur\R)jj^i + X 3 (sym(i?^ Curli?))^^^ 

+ pc{Rji,j ~ Rij,j)- (2-16) 
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We can also write this in matrix notation. Using the matrix divergence, we have 


p dttVL = (m + Me) Au + {p + \- Pc) grad div u 

+ xi grad tr(i?^ Curl R) + xs Div[dev sym(i?^ Curl R)] 

— 2/ic Div [skew i?]. (2-17) 

This is the equation of motion for the displacement u. The corresponding equations of motion for the 
rotations are considerably longer. 

Alternatively, we can write (12.171) using a slightly different notation involving the matrix Vu. We 
have the following identities 

p Au + (/r + A) grad div u = Div [2^sym Vu + Atr(Vu)l], (2.18) 

Pc Au — Pc grad div u = 2/Xc Div [skew Vu]. (2-19) 

Therefore we can write (12.171) as follows 

pdttu= Div [2/rsym Vu + Atr(Vu)l] 

+ xi grad tr(i?^ Curl R) + xs Div[dev sym(i?^ Curl R)] 

+ 2pc Div[skew Vu] — 2pc Div[skew i?]. (2.20) 

Lastly, using that skew(i? — 1) = skewi? we arrive at 

pdttu = Div [2/rsym Vu + Atr(Vu)l + 2/icSkew(Vu — {R— 1))] 

+ Xi grad tr(i?^ Curl i?) + xs Div[devsym(i?^ Curli?)], (2-21) 

which is a neat form of the equation for the displacements. 


3 Solitonic solutions 


3.1 Assumptions 

The full equations of motion for the rotations cannot be explicitly written down since they would be 
too long. We can simplify our equations considerably by making several assumptions, most of which 
are related to the ones used in [3]. 

First of all, let us assume that the points in our continuum can only experience rotations about 
one axis, the z-axis say, which means we choose 


/O -ct> 
A = {(j) 0 

yo 0 


o\ _ 

0 , R = 

0 / 


( cost/) 
sin^ 
0 


— sin (j) 0 \ 

cost/) 0 , 

0 l) 


(3.1) 


where (j) = (j){x^y, z,t) is an arbitrary function of the three spatial variables and time. Our choice of 
the rotation axis sets both oi and 02 to zero. Analogously, we assume that the displacements must 
occur along the axis of rotation. Therefore, we have that ui and U 2 vanish, and M 3 = y, z, t), with 
'0 also being an arbitrary function of the three spatial variables and time. 

It is well known that there are two possible types of plane waves in classical elasticity: transverse 
and longitudinal waves. Since we assumed that our wave propagates along the z-axis, we see that if 
the elastic wave is longitudinal, we have that ip = 0(z,t). On the other hand, if the elastic wave is 
transverse, then tp should only depend on the coordinates orthogonal to the z-axis, i.e. '0 = ''P{x,y,t)- 
While it is less clear what it means for a rotational wave to be transverse or longitudinal, we will simply 
use the above conventions. For a longitudinal wave we would write 0 = 0(z,t) while for a transverse 
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rotational wave we assume (j) = 4>{x-, y, t). Since we are interested in finding wave type solutions of the 
equations of motions, we can in principle study four different combinations of the above waves. 

In the following we will discuss the equations and their solutions assuming that both, the elastic 
and the rotational wave are both longitudinal waves, this means we assume t) and (j) = (j){z, t), 

or expressed in our original variables 


u = 



\0 0 l + dzipj 


for the desplacements. For the rotations we have 


(3.2) 


f cos (j){z,t) 

— sin^(2;, t) 

o\ 

rj. fdzCp 

0 

o\ 


sin(/)(z,t) 

cos (j){z,t) 

0 , 

R Gurli? = 0 

dzfp 

0 . 

(3.3) 

V 0 

0 

1 / 

V 0 

0 

0/ 



This ansatz for the solution can be interpreted in the following way. Along the z-axis we have a 
longitudinal wave which we can view as a compression wave. This axis also defines our axis for 
rotations which means we consider rotational waves transversal to the compression wave. 


3.2 The equations of motion 

Using our ansatz, the six independent equations of motion reduce to two equations, one for the field (f> 
and one for the field tp. They are given by 


-(ki + 6K3)d^^4>+ -(3xi - X3)dzztlJ - Me sm{(j)) - Protdtt4> = 0, 

O D 

(A + 2y)dz,^ip + -(3xi - X3)dzz(p - pdttpJ = 0. 

The two equations can be rewritten, using matrix notation, as follows: 




Me /sin(^)\ 

Prot \ ^ / 


where the matrix M is given by 


/ {ki + 6k3)/{3p 

rot) 

V 2(3X1 - X3)/(3p) 


(3X1 - X3)/{Gpvot)'\ 
{X + 2 p)/p J 


(3.4) 

(3.5) 


(3.6) 


(3.7) 


We note that the component M 22 is the square of the elastic (longitudinal) wave speed and we will use 
the notation M 22 = I'eias- Similarly, the component Mu corresponds to the square of the rotational 
wave speed and we will write Mu = It turns out to be convenient to work with the components 
of M instead of the elastic moduli of the energy functional. 

The system of equations (I3.4I) - (I3.5I) shows similarities with the sine-Gordon-d’Alembert system, see 
for instance [H, which has been previously studied in the context of micropolar elasticity and gives 
rise to soliton-like solutions [li [111301131] ■ The key difference between our system and those studied 
previously is the presence of interaction terms due to the interaction Lagrangian (12.31) which results in 
additional second derivative terms. Moreover, the sine-Gordon-d’Alembert system typically contains 
first derivative terms of the field variables which are also absent is our approach. It is likely that a 
more general ansatz will result in equations of motion which show a greater similarity to the sine- 
Gordon-d’Alembert system, however, the interaction terms due to interaction will make this system 
uniquely different again. 
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3.3 Sine-Gordon type equation 

We start by recalling the scalar wave equation x) — c^dxx‘f(t, cc) = 0 for the function ip(t^ x) where 

c is the wave speed. A very important and related equation is the sine-Gordon equation dtt‘^(t^x) — 
dxx^fit, x) + sin((/j(t, a;)) = 0 which is known to have soliton solutions. 

Inspection of (13.611 motivates us to seek soliton type solutions of this coupled system of equations. 
Equation (13.41) contains the term (Syi — xs)/{^Prot)dzz'4’- If we can ‘remove’ this term we would obtain 
a sine-Gordon equation for </>. The easiest way to obtain a sine-Gordon equation would be to simply 
set xi = X 3 = 0 in which case the displacement = ip will satisfy the wave equation, while the 
rotation pi will be of the form of a sine-Gordon equation. This rotational sine-Gordon equation was 
discussed in [2]. 

In the following we will demonstrate how this can be achieved. Start with a plane wave ansatz for 
Ip in the form 


ip = g{z- vit), 


(3.8) 


which corresponds to a travelling wave in the positive z direction only, and g is an arbitrary function 
of its argument, and vi is a wave speed. Due to the complicated structure of the equations, we were 
not able to consider the more general wave ansatz f{z -I- vit) -|- g{z — fit). Now that ip satisfies the 
wave equation dttip = v\ dzzi’, this is substituted into (13.51) to obtain 


g'\z - vit) = 


M- 


21 


-dzz(p, 


(3.9) 


^elas 


where we assume that vf — 022 0. Next, relation (13.91) is used so that dzz'P’ can be eliminated 

from (j3.4L we thus arrive at 


o j, 12, -^12 M21 
OttP - ( frot + “2 Z2 — 


^elas 


dzz(p H- —sincp = 0, 

Prot 


(3.10) 


which is the sine-Gordon equation. It could be brought into the above mentioned standard form by 
rescaling the spatial coordinate 2 ; as follows 


z = [ V, 


Mi2 M 21 \ 


''elas 




ft - f: 

and introducing the parameter rr? = /ic/prot, and we would arrive at 

dttp> - dii<p -I- rr? sin ^ = 0. 

However, in the following we will work with our original variables. 


(3.11) 


(3.12) 


3.4 Solitonic solutions 

We are now seeking solitonic solutions for the function (p and assume it be of the form 

(^ = 4arctane'=("-"=‘)+^ (3.13) 

where V 2 , k and S are constants. This sine-Gordon ansatz is substituted into Eqs. (13.91) with the aim 
of solving for g{z — fit). In order to proceed we must make one additional assumption, namely that 
the elastic and the rotational wave propagate with the same wave speed f 1 = f 2 = f. By substituting 
the ansatz (|3.13l) along with s = z — v t into (13.51) we arrive at the following differential equation for 
ais) 

4M2i/c2 eka+S(^^2iks+S) 

^ ^ {e^{ks+S) + 1)2 


(3.14) 









which can be integrated twice with respect to s to give 


g{s) = - 4 M 21 Ci+ C2S, (3.15) 

^elas - « 

where Ci and C 2 are two constants of integration which we set to zero. 

Next, we need to substitute the ansatz (13.131) with (13.151) into the remaining equation (13.41) . This 
yields an algebraic equation which relates the various constants of the problem. We find, after some 
algebra, 

— iv{M)v^ + det(M)) — (feias ~ = 0, (3.16) 

Prot 

which fixes either the parameter k or the velocity v, and we should note that tr(M) = The 

term det(M) also depends on the coupling constants. This relation for k is given by 




^elas - « 


— ti{M)v^ + det(M) Prot' 


(3.17) 


The solution with positive root is called a kink, while the solution with the negative root fc_ is 
known as the antikink, see 

We can also use equation (13.161) to find in terms of fc, m and each Mij. As expected, letting 
Pc = 0 in (13.161) gives us an equation which determines the velocities independent of k. The resulting 
velocities are those reported in equation (4.18) of [5], since this is equivalent to setting pc = 0. 
Putting everything together, our solution is given explicitly by 

()> = 4arctane'=±(^-’'‘)+^ (3.18) 

^ = -J^^HL^4arctane'=±("-"‘)+^ (3.19) 

^elas - ^ 

with k± given by (13.171) and <5 and v being constants. 



Figure 1: Visualisation of the soliton solution. The material points are represented by tori, the 
microrotations are visualised by the rotations of the tori. 

We visualise this solution in Fig. [TJ The z-axis is along the horizontal. Every material point with 
microstructure is represented by a ring (torus). The grey pattern on the rings helps to visualise the 
rotations about the z-axis. Looking closely at the figure, we see that rings 1 & 2, and 5 & 6 are in 
their (almost) undeformed states while the positions of rings 3 & 4 are displaced. Firstly, we note 
that the spatial separation between the middle rings is smaller than the one of the outer rings. This 
displacement is due to the longitudinal elastic soliton travelling along that axis. Secondly, we note that 
the middle two rings are rotated with respect to the outer rings. This visualises the microrotational 
soliton of the microstructure. 
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4 Discussion 


The geometrically nonlinear Cosserat micropolar elasticity model with nonstandard constitutive cou¬ 
pling between Cosserat rotations and the deformation gradient was studied. We were able to find 
soliton type solutions assuming small displacement but allowing arbitrarily large Cosserat rotations. 
In general the Euler-Lagrange equations of this system are very complicated, however, using a specihc 
ansatz, it is possible to significantly simplify these equations. We assume a compression wave along 
the z-axis and restrict possible rotations to be about this axis. The corresponding solution we find 
consists of two solitonic waves, an elastic one and the microrotational one, both of which have the 
same wave speed. 

The existence of such soliton type solutions in our model is related to the Cosserat couple term. In 
the limit 0 our equations of motion would not have such a solution, and we would instead find two 

coupled plane waves. It is precisely this Cosserat couple modulus which appears to localise these 
waves. At this stage it is not clear whether this statement can be generalised to wave type solutions 
of different form, our ansatz is very specific. However, these results point into a possible direction of 
interpreting the Cosserat couple modulus. 

The dynamical Cosserat model studied here is particularly interesting from various points of view. 
The model is intrinsically nonlinear because rotations in three dimensions do not commute, in other 
words the group SO (3) is non-Abelian. In the case of planar Cosserat elasticity, the matrix rotation 
R has one degree of freedom only, one can only rotate about one axis, and these rotations commute. 
The group SO (2) is Abelian. Hence, one cannot meaningfully simplify the dynamical Cosserat model 
to the planar case as essential features of the model would be lost. It will be important to study the 
complete model further, clearly the next significant task is to study the fully nonlinear model where the 
displacements are not assumed to be small. This means studying a continuum with arbitrarily large 
displacements and arbitrarily large rotations. One should think that a suitable ansatz will simplify 
these equations to the point where one can seek explicit solutions. 
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